*------------------------------------------------------------------------------*
					*** BOXPLOT PM10 BY DECILES ***
*------------------------------------------------------------------------------*

use "${tempdata}kreise_complete_all.dta", clear

keep if samplePM10_allcontrols == 1

*------------------------------------------------------------------------------*

***Divide counties into deciles based on average levels of PM10
egen AV_PM10_10ug = mean(PM10_30km_10ug), by(kreis_code)
xtile deciles_AV_PM10 = AV_PM10_10ug, nq(10)

*------------------------------------------------------------------------------*

***Box-Plot
graph box PM10_30km_10ug, over(deciles_AV_PM10) ytitle("PM10 (10 {&mu}g/m{sup:3})") b1title("Average level of PM10 (decile)") graphregion(color(white))

graph export "${outputs}figure_S1b.png", replace

*------------------------------------------------------------------------------*

clear

exit
